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We study the effect of a magnetic field on the current-phase relation of a topological Josephson 
junction formed by connecting two superconductors through the helical edge states of a quantum 
spin-Hall insulator. We predict that the Zeeman effect along the spin quantization axis of the helical 
edges results in an anomalous Josephson relation that allows for a supercurrent to flow in the absence 
of superconducting phase bias. We relate the associated field-tunable phase shift 0o in the Josephson 
relation of such a 0o-j unction to the existence of a so-called helical superconductivity, which may 
result from the interplay of the Zeeman effect and spin-orbit coupling. We analyze the dependence 
of the magneto-supercurrent on the junction length and discuss its observability in suitably designed 
hybrid structures subject to an in-plane magnetic field. 
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I. INTRODUCTION 


The topological properties of quantum spin-Hall in¬ 
sulators (QSHI) manifest themselves in current-carrying 
helical edge states, characterized by a locking of their 
group velocity to the spin orientation 1 The effective 
one-dimensional (ID) superconductivity induced in those 
edge states by conventional superconductors (S) via prox¬ 
imity effect is predicted to also be topologicaP. This 
is expected to lead to a fractional Josephson effect in 
S-QSHI-S junctions, due to topologically protected An¬ 
dreev bound states (ABSJP 1 . Recent experiments on 
I IgTe/CdTcP or InAs/GaSb 1 - as QSHI, contacted with 
conventional superconducting leads, demonstrated that 
the Josephson current is mainly carried by edge states in 
the QSH regime. However, a clear signature of the topo¬ 
logical superconductivity induced on edge states is still 
lacking. 

The role of a Zeeman field transverse to the spin quan¬ 
tization axis of the edge states has been discussed in¬ 
tensively. In the “bulk” it may induce a transition from 
topological to topologically-trivial superconductivitjPH 
whereas a local field in the junction area acts as an ef¬ 
fective barrier®. Much less is known about the role of a 
Zeeman field parallel to the spin quantization axis. An 
anomalous Josephson effect, that is, a supercurrent flow 
when no superconducting phase bias is applied to the 
junction, was calculated in S-QSHI-S junctions with a lo¬ 
cal Zeeman field applied in the junction aresPP Similar 
results in ju nction s with 3D topological insulators were 
also obtaine d 13 * 14 !. These systems realize so-called 00- 
junctions, where the current-phase relation has a phase 
shift 0o that is tunable with the external magnetic held. 
The effect found in those works disa ppea rs in the short- 
junction limit. However, experiments 9 10 involve a mag¬ 
netic Held applied to the entire system. We show that, in 
this case, the anomalous Josephson effect depends both 
on the Held in the superconductor and the field in the 
junction area via two different though related mecha¬ 


nisms. As a result, we find an anomalous current in 
junctions of any length, with an amplitude that can be 
as large as their critical current in the absence of a mag¬ 
netic field. Specifically, we show that the effect of the 
field in the superconductor allows one to probe whether 
superconductivity is indeed induced in the edge states. 

Note that an anomalous Josephson effect may occur 
in any junction with broken time-reversal symmetrjG®. 
Indeed, 0o-junctions were discussed previously in a va¬ 
riety of topologically-trivial systems, mainly in the pres¬ 
ence of spin-orbit coupling and a Zeeman field applied 
to the junction are£p®®3, bu t als o in superconductor- 
ferromagnet hybrid structures 23 24 . In the former case, 
the predicted effect is small. Namely, when both spin- 
orbit-induced helical bands cross the Fermi level, a par¬ 
tial compensation between the two helicities occurs and 
only a residual effect proportional to the mismatch in 
their densities of state remains. In particular, in a 
ID system, such a mismatch requires deviations from a 
parabolic spectrum^. The predicted anomalous Joseph¬ 
son effect has not been observed experimentally in these 
systems, so far. The large anomalous current we find in 
S-QSHI-S junctions is a direct consequence of the heli¬ 
cal nature of the edge states. It should be observable in 
setups similar to those studied in Refs. [9l and ilOl 


The outline of the article is as follows: We introduce 
the model in Sec. [TTJ Then, we study the anomalous 
Josephson effect in a S-QSHI-S junction along a single 
edge in Sec. Ill In particular, we determine the current- 
phase relation as a function of the external magnetic field 
and the junction length. We turn to the observability of 
the effect in S-QSHI-S junctions where both edges con¬ 
tribute in Sec. [iVj We argue that Josephson junctions of 
unequal lengths should be realized in order to observe the 
effect. Finally, in Sec.|Vj we discuss that the anisotropy 
of the gyromagnetic tensor should allow for the observa¬ 
tion of the effect with an in-plane magnetic field. We also 
show that the effect is robust with respect to a small mis- 
alignement between the applied field and the spin quan- 
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tization axis. 


II. THE MODEL 

A Josephson junction formed along one of the edges 
of a QSHI can be desc ribed by the ID Bogoliubov-de 
Gennes HamihoniarPG 1 


H = ( v F p x a 3 - p)r 3 - ha 3 + A(x)r + + A* (a;)r_ . (1) 


Here, vp is the Fermi velocity, x and p x are the position 
and momentum operators, respectively, /i is the Fermi en¬ 
ergy, and h > 0 is a Zeeman field along the spin quantiza¬ 
tion axis. (The effect of a perpendicular component of the 
field will be discussed below.) The superconducting gap 
induced by conventional superconducting leads is given 
as A(x) = A o [e"^/ 2 0(x - L/2) + e^/ 2 0(-x - L/2)], 
where A 0 is the magnitude of the gap, <f> is the phase 
difference between the two leads, and L is the junction 
length. Moreover, cq and r* are Pauli matrices act¬ 
ing in spin and Nambu spaces, respectively, and t± = 
(ri ± 1 T 2 ) /2. Here we use units where h = kp = 1. 

The role of the magnetic field within the supercon¬ 
ducting regions is most easily understood by consider¬ 
ing first a ID “bulk” superconductor, i.e., by setting 
L = 0 and <j> = 0 in Eq. (jTJ) . The Zeeman term in¬ 
duces a momentum mismatch 2 h/vF between left- and 
right-moving states at the Fermi level which may be 
gauged out using the unitary transformation H -A U^HU 
with U = e *(VvF)zT 3 . However, this gauge transforma¬ 
tion modifies the order parameter, Ao —> Aoe~ 2l ^ h ^ VF ^ x . 
Thus, for a uniform order parameter, A 0 = const, one 
obtains a current-carrying excited state, whereas the 
ground state would require a spatially modulated order 
parameter with wavevector q, Aq(x) = Aoe tgx . 

Indeed, the free energy density of the system can be 
easily computed and depends on the modulation wavec- 
tor q through an effective field h q = h — VFq/%- For 
details, see the appendix. At zero temperature, one finds 
F(h,q ) = F 0 + A%/ (2ttvf) f {h q /A 0 ), where F 0 is inde¬ 
pendent of q and /(x) = x 2 + 0(|x| — l)[arccoshx — 
\x\\/x 2 — 1]. The supercurrent is obtained using the ther¬ 
modynamic relation / = —2 e(dF/dq). One readily shows 
that the free energy density is minimized and the cur¬ 
rent is zero for h q = 0, corresponding to a modulation 
wavevector q = 2 h/vp- Such a modulated or so-called 
“helical” superconductivity has been studied in higher 
dimensions^. 

By contrast, if superconductivity is induced by a con¬ 
ventional bulk superconductor with constant phase, the 
induced order parameter inherits the bulk superconduct¬ 
ing phase, and a modulation is not possible. Then q = 0, 
and the superconductivity induced in the edge states car¬ 
ries a current, 


m = - 

7 r 


h-Q(h- Ao)^ 2 - A 2 


( 2 ) 


This is precisely the current /(</> = 0) that would flow in 
a short junction, L <C £ with £ = vf/ Ao, at zero phase 
difference. Thus, the fact that the proximity-induced su¬ 
perconductivity forces the system into an excited state 
yields an anomalous Josephson effect. The anomalous 
current increases proportionally to h at h < Aq and then 
decreases as / ~ eAg/(27 ih) at h A 0 . 

In the following, we extend the result ([ 2 ]) to arbi¬ 
trary junction lengths and temperatures, and study the 
current-phase relation. Note that the fact that Ao is an 
induced gap also implies that there is no self-consistency 
condition and that fields h > Ao are possible as long as 
Ao is sufficiently smaller than the intrinsic gap of the 
superconducting leads. 


III. THE CURRENT-PHASE RELATION 


We use the formalism of Refs.^^H to obtain the 
Josephson current from the Hamiltonian 0, 

1 - a 2 (uj u - ih)e- 2 ^- ih y EL e^ 

i /—0 

( 3 ) 

Here = (2v + l)nT are Matsubara frequencies at tem¬ 
perature T, a(w) = i(uj — \Jw 2 + Ag)/A 0 , and Ep = 
vf/L is the Thouless energy of the junction. Equa¬ 
tion © accounts for the contributions of both the states 
in the continuum outside the superconducting gap and 
the ABS, whose subgap energies E n correspond to the 
poles of the r.h.s of Eq. ^ after analytic continuation, 
—» — iE + 0 + . In particular, the ABS energies read 


Ij = -4eT—SRy In 
dm ^' 


2 arccos 


/ E n + h \ 

V^T J 


2 (E n + h) 
E~l 


(f> + 27m, n £ Z. 


Equation ([ 3 ]) can be used to numerically compute the 
current-phase relation at arbitrary junction lengths and 
temperatures. The results at low temperatures and vari¬ 
ous fields are shown in Fig. [I] The current-phase relation 
and the corresponding anomalous Josephson current as 
a function of the magnetic field are shown for a short 
junction [L = 0.1£, panels (a)-(b)] and a long junction 
[L = 10£, panels (c)-(d)], respectively. Below we analyze 
both short and long junctions further, starting with the 
limit of zero temperature. 


A. Short vs long junction 

In the short junction limit, Ao Ep, we find that the 
continuum states are essential in determining the current- 
phase relation (in contrast with conventional short junc¬ 
tions, where the supercurrent is carried by ABS onljl 27 ^). 
Evaluating Eq. ^ at </> = 0, one readily recovers the 
result © which is a pure continuum contribution. At 
finite q i, the junction accommodates for a single bound 
state with energy Ea = A o cos(0/2) — h. The unique 
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zero-energy solution at <j>* = 2 arccos(/i/A 0 ) for h < A 0 
is a consequence of the topological nature of the junc¬ 
tion. It leads to a jump in the current phase relation, 
cf. Fig. 0a), which disappears at h > A 0 , signaling the 
transition to a topologically trivial state. The current- 
phase relation can be obtained by expanding Eq. ([3]) in 
harmonics and evaluating each term. Summing up the 
series, one finds, at T = 0 and h < A 0 , 


h ) 


eh eAo . (f> . 

-b „ sin - sign 

7T Z Z 


sin 



( 5 ) 


where the two terms correspond to the continuum and 
the ABS contributions, respectively®]. The bulk contri¬ 
bution due to the field in the superconductors yields an 
asymmetry between the critical currents in opposite di¬ 
rections, /+>/“. Such an anomalous Josephson effect 
is, thus, a direct probe of the nature of the induced su¬ 
perconductivity underneath the contact^! 

We now turn to the long junction limit, Ao El- At 
h < A 0 , Eq. (|4f yields a large number of ABS with en¬ 
ergies E n = —E l [(/)+ (/) h + 2n(n + 1/2)] /2, where fa = 
2/i/£x +2 arcsin(/i/Ao). The phase shift fa has two con¬ 
tributions: The first term is proportional to the junction 
length and can be traced back to the magnetic field in 
the junction area, it is dominant for El <C A 0 . The 
second term stems from the bulk effect discussed above. 
Correspondingly, the current-phase relation at T = 0 is 




eE L 

2n 


<j> + fa- 27TInt 


( 6 ) 


which extends the result obtained for lon g S-Q SHI-S 
junctions in the absence of a magnetic fiel d 32 * 33 ! The 
anomalous Josephson current, thus, displays a sawtooth 
behavior as a function of the applied magnetic field, 
which is visible in Fig. 0d). This is reminiscent of the 
Little-Parks effectP^, though with a paramagnetic rather 
than orbital origin. As in the short junction limit, the 
topological nature of the junction manifests itself in a 
jump in the current-phase relation when the lowest ABS 
reaches zero energy, at <j>* = (n — fa) mod 27r. 

At larger fields, h > A 0 , one finds 



FIG. 1. Anomalous Josephson effect in short [L = 0.1£, 
panels (a) and (b)] and long [L = 10£, panels (c) and (d)] S- 
QSHI-S junctions. Panels (a) and (c) show the current-phase 
relation at temperature T/Aq = 10 -3 for different values of 
the applied magnetic field h. Panels (b) and (d) show the 
anomalous Josephson current at 4> = 0 as a function of h. 


function of the field near h = 0 is the same in junctions 
of any length. 

Namely, 


= 0 ,h) 
dh 


h=0 


e 
7r 


( 8 ) 


This universal result is a consequence of the helical nature 
of the QSHI edge states. 


B. Finite temperature effects 


7/(0, h ) = 


eE L 
7r 


arc tan 


sin (</>+f/) 


2 arccosh -^ L - 

e A o 


COS {(j) + 



( 7 ) 

As expected, there is no more jump in the current-phase 
relation, and the anomalous Josephson current is sup¬ 
pressed with increasing field. 

The behavior of short and long junctions is quite differ¬ 
ent. In short junctions, the anomalous Josephson effect 
stems from the magnetic fields in the leads. By con¬ 
trast, in long junctions, the dominant contribution at 
small fields comes from the magnetic field in the junction 
area. In general, both the field in the leads and in the 
normal part of the junction contribute. Note, however, 
that the slope of the anomalous Josephson current as a 


Finite temperatures smear out the sharp features in the 
current-phase relation. At T> min[Ao, El], a sinusoidal 
behavior, 

Ij(4>,h) = I C (A 0 , h,T, E L ) siny> + fa(h,T,E L )], (9) 

is found for both short and long junctions. In particular, 
for short junctions, we find I c = eAg |V ; i(“)| /(47r 2 T) and 
fa = arg{'0i(z)}, where fa is the digamma function and 
z = 1/2 — ih/(2nT), such that the phase shift increases 
from 0 to 7 t/2 with the field. In long junctions, we find 
I c = 4 eT exp[—27rT/£x]|a(27rT— ih)\ 2 and fa = 2 h/EL + 
2arg{a(27r T — ih)} + n. For h <C T <C Ao, one obtains 
\a{2-KT—ih)\ = 1 and arg{a(2-7rT — ih)} = — 7r/2, whereas 
for T <C A 0 <C h, one obtains |a(27rT— ih)\ = Ag/(2 h) 2 
and arg{a(27rT — ih)} = 0. 
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IV. DOUBLE JUNCTIONS 

When creating a Josephson junction with a QSHI, typ¬ 
ically both edges of the QSHI contribute to the Josephson 
currenlPESl. If the width of the QSHI is sufficiently large, 
the system may be described as two junctions in paral¬ 
lel and their contributions may be computed separately. 
Then, as the two edges of the QSHI have opposite helici- 
ties, the contribution of the second edge can be accounted 
for by another copy of Hamiltonian (IH with h —> —h. 
The corresponding current is Ij(4>,—h) = — Ij(—cf>,h ). 
Similar to the case of conventional 00 -junctions discussed 
in the Introduction, adding the current contributions 
of the two edges leads to a (partial) compensation of 
their anomalous Josephson currents. Here, however, the 
spatial separation of the two helicities makes an im¬ 
portant difference. Only if the two junctions on either 
side of the sample have the same length the compen¬ 
sation is exact, and we obtain the conventional result, 
jt ° tal (0 = Q fo'j — o However, if the two junctions have 
unequal lengths, as shown schematically in Fig. |2j)a) , the 
compensation is only partial and a residual effect re¬ 
mains. This residual effect is a signature that the Joseph¬ 
son current is carried by helical edge states. The depen¬ 
dence of the anomalous Josephson current at 0 = 0 is 
plotted as a function of h for various temperatures in 
Fig. §b). Note that is straightforward to take into ac¬ 
count the additional phase shift between the two edges 
due to the orbital effect of the field, if the junction area 
encloses a magnetic flmf 9 10 35 l 

V. DISCUSSION AND CONCLUSIONS 

We now turn to the conditions of applicability of our 
model 0. According to the Bernevig-Hughes-Zhang 
(BHZ) model for inverted electron-hole bilayeriP, the spin 
quantization axis points along the growth direction of the 
heterostructure, and the Zeeman field needed to obtain 
a 0 o-junction would originate from an out-of-plane mag¬ 
netic field. This configuration would most likely suppress 
superconductivity in the leads. However, in real systems, 
the BHZ model should be supplemented with bulk inver¬ 
sion asymmetry (BIA) and Rashba spin-orbit coupling 
terms™ Those terms result in a tilt of the quantization 
axis toward the quantum well plane, thereby allowing for 
the generation of the Zeeman field appearing in Eq. 0 
with an in-plane magnetic field. 

To obtain the dependence of this Zeeman field on an 
external magnetic field B, one needs to determine the 
effective Zeeman term in the Hamiltonian for the helical 
edge states, H z = )Tb =1 2 3 (t, • B)cr.j. Here Q are vectors 
that can be computed within the extended BHZ modeP2 
If the field is applied along the direction ti x t 2 , no spin 
gap opens in the edge excitation spectrum, and one ob¬ 
tains Eq. (|T|) with an exchange field h = — 13 • B. Spe¬ 
cific values for t 12 in a 7nm-thick HgTe/CdTe quantum 
well 3 ^ indicate that the preferential direction lies close to 



FIG. 2. (a) Proposed setup to detect the 0o-junction be¬ 

havior in a S-QSHI-S hybrid system: The magnetic held B 
is applied in the plane of the quantum well. The edge states 
on both sides of the sample contribute to the Josephson cur¬ 
rent. A net anomalous Josephson effect remains, if the junc¬ 
tions have unequal lengths, L± ^ L 2 . (b) The anomalous 

Josephson current at 0o as a function of the Zeeman splitting 
h = /rsg e fi|B|/2 in the case L 2 = 10£ and L\ = L 2 /3 for 
various temperatures. 



FIG. 3. The slope of the anomalous Josephson current in 
short junctions as a function of the magnetix held h± perpen¬ 
dicular to the spin quantization axis and the chemical poten¬ 
tial. Th e slope re mains close to its quantized value as long as 
h± < y/Ag +/r 2 . 


the plane and perpendicular to the edges. The effective 
5 -factor is expected to be fairly larg^S^. 

Experimental studies of the magneto-conductance 
anisotropy show indeed that the conductance is hardly 
affected by an in-plane field, while a perpendicular mag¬ 
netic field leads to a large suppression^. These re¬ 
sults also show that the topological protection against 
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backscattering, although in principle not guaranteed 
when time-reversal symmetry is broken by an in-plane 
magnetic field, is in practice approximately conserved. 

Our results are robust with respect to the opening of 
a small spin gap due to a misalignment between the field 
and the spin quantization axis. The anomalous Joseph- 
son current in a short junction for an arbitrary direction 
of the applied magnetic field is derived in the Appendix. 
In Fig. [3} we show dhlj(<j) = 0,/i)|h=o as a functions of 
the field h± perpendicular to the spin quantization axis 
and the chemical potential [i. The slope remains close 
to its quantized value as long as h± <C \/Ag + /r 2 . By 
contrast, when h± > y / A§”+/?, the system becomes 
topologically trivial, and the effect disappears rapidly. 

In conclusion, we have demonstrated that the helical 
nature of the QSHI edge states leads to an anomalous 
Josephson effect in S-QSHI-S junctions subject to a mag¬ 
netic field. The resulting anomalous supercurrent, flow¬ 
ing at zero phase difference between the two supercon¬ 
ducting leads, is field tunable. Both the field in the su¬ 
perconductor and in the junction area contribute to the 
effect and probe the helical nature of the edge states in 
the corresponding parts of the system. We also discussed 
how to observe this effect using hybrid structures based 
on available QSHI realizations, analyzing the contribu¬ 
tions of both edges, the required magnetic field direction 
as well as the stability of the effect with respect to a finite 
chemical potential and a misalignment of the magnetic 
field and the spin quantization axis. Similarly, we expect 
a pronounced anomalous Josephson effect in jun ctions 
based on nanowires with strong spin-orbit couplin ^ 38 l' 39 l 
when they are in the topological regime. 
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Appendix: Free energy and anomalous Josephson 
current of short junctions 

In this appendix, we derive an expression for the 
anomalous current, 1(h) = Ij((f> = 0 , h), carried by one 
of the edges of a short S-QSHI-S junction for an arbitrary 
orientation of the Zeeman field, h = he 3 + h±e 3 , where 
ei and e 3 are perpendicular unit vectors, and the spin 
quantization of the helical edge state is along e 3 . 

As we argued in the main text, the anomalous cur¬ 
rent is obtained as I = —2e(dF/dq)\ q= o, where F is the 
free energy density of a “bulk” ID superconductor with 
a spatially modulated superconducting order parameter, 


A(a;) = Aoe* 9 *, where q is the helical modulation wavec- 
tor. In the following, we assume q > 0 for definiteness. 

As the linear spectrum of the helical edge states is not 
bounded from below, the free energy will depend on the 
large momentum cut-off, k c . To properly treat this cut¬ 
off, we write the Hamiltonian in the form 

T-L = "^2 \i v Fk - n - h)a\a k - (v F k + H - h)b\.b]lA.l) 

\k\<k c 

—h±(a\.bk + h\ak) 

I fc| <fce— f 




A 0 a 


k-\- 9 -k-\-T) 


+ h.c. . 


Here, a k and b k are annihilation operators for right- 
moving (spin up) and left-moving (spinjdown) electrons 
with momentum k and Fermi velocity vM Without loss 
of generality, we choose A 0 to be real and positive. As 
the pairing term couples states with different momenta, 
the effective cut-off becomes ^-dependent. Thus, when 
introducing the Bogoliubov-de Gennes Hamiltonian (cf. 
Eq. (1) in the main text), additional cut-off dependent 
terms have to be kept. Namely, the Hamiltonian (A.l) 
may be split into two parts 


n = 


^ rlH k r k + n > , 


2 ^ 

ifci<fc, 


(A.2) 


where 


H k = (v F ka 3 - h)t 3 + A 0 n - h±a 3 - h q a 3 , (A.3) 


with h q = h—v F q/2, is the Bogoliubov-de Gennes Hamil¬ 
tonian, which is expressed with the help of Pauli matrices 
< 7 j and n acting in spin and Nambu spaces, respectively, 

and 0 — (afc+^.G+j.bfc+^.-A+j) - Furthermore, 


■*> = -]»> XI <4 4) WfeV 4 ) 

td—± vk c <k<i>k c +q/2 ^ 

2 ^ 

|fc|<fc c 

with 

*^(■'‘7-*-,,^^) (*■») 

is the contribution from energies close to the cut-off, 
where modifications of the spectrum due to supercon¬ 
ductivity as well as the transverse field are negligible. 
Its contribution to the free energy density is obtained by 
taking the expectation value of 'H > in the ground state 
of the system, yielding 


F> 


Y ( v F k — fi — h) 

-k c <.k<—k c -\-q/2 

^ y (~v F k - n + h) 

k c <k<k c -\-q /2 
1 (/^ _ 


2ttv f 


(A. 6 ) 


yk c 

7r 
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Let us now turn to the Bogoliubov-de Gennes Hamil¬ 
tonian H k . It has four eigenenergies, Ei(k) with i = 
1, 2, 3,4, determined through the equation 

(. E 2 - \ 2 +k )(E 2 - X 2 _ k ) + 8 iiv F kh q E = 0 (A.7) 

with 


\± k = vpk 2 + n 2 + h 2 + h 2 ± + (A.8) 

±2yJv* F k*(f + hi) + ( M 2 + A g)(ft2 + hi). 

The general expression for the free energy density reads 
^=(l/i)E| fc |< fcc E i;Bi >o^[/T(^)-l/2]+F>, where 
fr(E) is the Fermi function at temperature T. Note that 
the solutions obey the relation Ei(k) = —Ei(—k). Thus, 
at zero temperature, the free energy density of the system 
takes the form 

k c 

F = -j n J dk J2 m k )\+ F> - ( a . 9 ) 

o 1 

The current is then given as 

k c 

1(h) = J dk J2dh\E,m + e -h, (A.10) 
o * 


where Ei(k) are the eigenenergies of the system at q = 0. 
Analyzing Eqs. (A.71 and(A.8), we note that zero-energy 
solutions exist if h z + h \ = Ag + p, 2 or \h\ > Aq. The 
topologically non-trivial region corresponds to low fields, 
h 2 + h\ < A 2 + /i 2 and \h\ < A 0 . 

We further analyze the result (A.101 in two limiting 
cases, namely at h± = 0 and for arbitrary h± in the 
limit h —> 0. 


1. Supercurrent for a Zeeman field parallel to the 
spin quantization axis 


When the Zeeman field is parallel to the quantization 
axis, h± = 0, the eigenenergies are 

E Sl s 2 (fe) = S 2 \J\v F k + Si/x ) 2 + Ag - Si h q , (A.11) 

where si, S 2 = ±1. Evaluation of Eq. (A.9) at zero tem¬ 
perature then yields 


F = F 0 + 


1 


27 TVf 


(A.12) 


+0(1 hq\ — Ag) 


Ag arccosh - \h q \^ h 2 q - A 2 


with 


Fn = -: 


1 


2-KVp 

+A 


(vf k c + / r ) 2 


(A.13) 


i+in 


/ 2 vpk c 

V^T 


that does not depend on q , in agreement with the expres¬ 
sion given in the main text. Evaluating the anomalous 
current 1(h) using Eq. (A.10), we then obtain Eq. (2) in 
the main text. In the topological regime, h < Ag, the 
large anomalous current is given as 1(h) = (e/n )h. 


2. Supercurrent response to a small Zeeman field 
along the spin quantization axis 

When the field component along the spin quantiza¬ 
tion is small, h 0, we may evaluate the eigenen¬ 
ergies perturbatively in h q . At h q = 0, we obtain 
Es°l 2 (k) = s 2 E Sl (k) with 


E Sl (k) 


^Upk 2 + fj, 2 + Ag + h\ 


(A.14) 



Making use of perturbation theory to obtain the correc¬ 
tion 5E SlS2 (k) to the energy E^\ 2 (k), we obtain the free 
energy density up to quadratic order in h q , 


F = F 1 


2ttvf 


(A.15) 


- vph\ / dk 


—k c 


[E + (k) + E.(k)] 2 ~^ 2 + Al) 
E+(k)E_(k)[E + (k) + E-(k )] 3 


where Fi does not depend on q. As the integral converges 
at large momenta, we may take the limit k c —¥ oo. The 
slope of the anomalous current in the limit h — > 0 is, 
thus, 


dh 


h—0 



—vfE 


(A.16) 


(E + (k) + E_(k)} 2 ^A(^ + A 2 ) 
E+(k)E_(k)[E+(k) + E_(k )} 3 


The dependence of the slope (A.16) on the chemical po¬ 
tential /i and the transverse field h± is shown in Fig. 3 
of the main text. A large signal, (dI/dh)\h=o ~ e/n, 
is obtained deep in the topological region, when h\ -C 

A§ + ft 2 - 
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